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Abstract 

We give analytic approximations to the baryon asymmetry produced by ther- 
mal leptogenesis with hierarchical right-handed neutrinos. Our calculation in- 
cludes flavour-dependent washout processes and CP violation in scattering, and 
neglects gauge interactions and finite temperature corrections. Our approxi- 
mate formulae depend upon the three CP asymmetries in the individual lepton 
fiavours as well as on three fiavour-dependent efficiency factors. We show that 
the commonly used expressions for the lepton asymmetry, which depend on the 
total CP asymmetry and one single efficiency factor, may fail to reproduce the 
correct lepton asymmetry in a number of cases. We illustrate the importance 
of using the flavour- dependent formulae in the context of a two right-handed 
neutrino model. 



Note Added: In the published version of this paper, we apphed incorrectly the A 
matrix which relates the asymmetries in 5/3 — to the asymmetries Y^a carried by 
the lepton doublets. This over-estimated the baryon asymmetry by a factor of order one 
(or a few). In this revised version, we modify sections 4.1.4 and 4.2 to give a more correct 
relation between our analytic estimates for the flavour asymmetries, and the final baryon 
asymmetry. We thank E.J. Chun and M. Plumacher for discussions on this issue, and in 
particular A Strumia for details of his calculations and results. 

Second Note Added: There are various other errors (we thank S Antusch, P Di Bari, 
EJ Chun, and C Gonzalez-Garcia for comments) which we correct in this version: 

1. In the previous net version (modified to use more correctly the A matrix), the overall 
sign of the baryon asymmetry was wrong in eqns and fl57p . 

2. The coefficient of the washout term in the Boltzmann Equations was a factor of 2 
too large: in the last term of eqns (ITT]) . ([IHD, and (IHUI) . and in the corresponding 
coefficients of the analytic solutions, the numerical factor should be 1/4 (not the 1/2 
of the published version) . 

3. The vev of the Higgs should be taken 175 GeV, not 246 GeV(given after eqn ([T])). 
This changes m* from the 2 — 3 x 10"'^ eV used in the published version (at eqns 
fl5|l . fl20|) ) to m^, ~ 10"'^ eV. NB: in the analytic approximations this erroneous factor 
of 2 cancels against the one of point 2 above. So the analytic approximations remain 
acceptable , using ^ 10^^ eV. 

4. Factors of two in eqns (1^21) . and ([2H]) have been corrected. 

5. In the published version, the captions of Figures 1,2, and 3 indicate that e was taken 
distinguishable from /i. However, for Mi > 10^ GeV, only the r is distinguishable. 
We therefore interchanged the e and r indices on the inputs {e^e — * e^r, 
Keey)- We also changed a washout factor in the caption of Figure 2 to account for 
points 2 and 3 above. We have remade the plots, with these modified inputs, and the 
plots remain unchanged. 

1 Introduction 

There are robust observational evidences that a tiny excess of matter over antimatter 
was produced in our Universe [1], but its origin is still a mystery. Baryogenesis through 
Leptogenesis [2] is a simple mechanism to explain this baryon asymmetry of the Universe. 
A lepton asymmetry is dynamically generated and then converted into a baryon asymmetry 
due to (i? + L)-violating sphaleron interactions [3] which exist in the Standard Model (SM). 

A simple model in which this mechanism can be implemented is "Seesaw" (type I) [Ij, 
consisting of the Standard Model (SM) plus two or three right-handed (RH) Majorana 
neutrinos. In this simple extension of the SM, the usual scenario that is explored (referred 

^unless tan /3 > 5 in a SUSY model 



2 



to as "thermal leptogenesis" ) consists of a hierarchical spectrum for the RH neutrinos, 
such that the lightest of the RH neutrinos is produced by thermal scattering after inflation, 
and subsequently decays out-of-equilibrium in a lepton number and CP-violating way, thus 
satisfying Sakharov's constraints. 

In recent years, a lot of work [5], El [7], has been devoted to a thorough analysis of 
this model, giving limited attention to the issue of lepton flavour [8]. The dynamics of 
leptogenesis is usually addressed within the 'one-flavour' approximation, where Boltzmann 
equations are written for the abundance of the lightest RH neutrino, responsible for the out 
of equilibrium and CP asymmetric decays, and for the total lepton asymmetry. However, 
this 'one-flavour' approximation is rigorously correct only when the interactions mediated 
by charged lepton Yukawa couplings are out of equilibrium. 

In ref. [H], flavoured Boltzmann Equations were written down. Flavour effects in "reso- 
nant leptogenesis" were studied in [9], discussed for thermal leptogenesis in the two-right- 
handed neutrino model in [TU], and used in [TT] to protect an asymmetry made in the 
decay of the middle right-handed neutrino. In the four generation models of [12], flavour 
was used to enhance the asymmetry, foreshadowing the results we obtain here. The impact 
of flavour in thermal leptogenesis has been recently studied in some detail [l3l[Tlj, including 
the quantum oscillations/correlations of the asymmetries in lepton flavour space [I3]. It 
was shown that the Boltzmann equations describing the asymmetries in flavour space have 
additional terms which can significantly affect the result for the final baryon asymmetry. 
In [13], we focused on how flavour effects can enlarge the area of parameter space where 
leptogenesis can work: the lower bound on the reheat temperature is mildly decreased, and 
the upper bound on the light neutrino mass scale no longer holdfl 

Flavour effects have not usually been included in leptogenesis calculations. This is 
perhaps because perturbatively, they seem to be a small correction. For instance, if the 
asymmetry is a consequence of the very-out-of-equilibrium decay of an initial population 
of right-handed neutrinos, then the total lepton asymmetry is of order e/(?*, where e is the 
total CP asymmetry in the decay, and counts for the entropy dilution factor. Clearly the 
small charged lepton Yukawa couplings have no effect on e. However, realistic leptogenesis 
is a drawn-out dynamical process, involving the production and destruction of right-handed 
neutrinos, and of a lepton asymmetry that is distributed among distinguishable flavours. 
The processes which wash out lepton number are flavour dependent, e.g the inverse decays 
from electrons can destroy the lepton asymmetry carried by, and only by, the electrons. 
The asymmetries in each flavour are therefore washed out differently, and will appear with 
different weights in the final formula for the baryon asymmetry. This is physically inequiva- 
lent to the treatment of washout in the one-flavour approximation, where indistinguishable 
leptons propagate between decays and inverse decays, so inverse decays from all flavours 
are taken to wash out asymmetries in any flavouiH. 



^The bound is removed when flavour efi^ects are relevant, which is the case for leptogenesis at tempera- 
tures < 1012 GeV. 

■^The "one-flavour" formulae describe leptogenesis that takes place at temperatures larger than 10^^ 
GeV, before the charged lepton Yukawas come into equilibrium. They are also appropriate for right- 
handed neutrinos who decay only to one flavour. (But note from eqns (|51l52p . that flavour effects can be 
important when there are small branching ratios to other flavours.) 
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In this paper we provide the necessary analytical expressions for the computation of the 
baryon asymmetry including flavour that the interested reader may apply to their preferred 
model. By comparing to the usually adopted 'one-flavour' approximation, we will show that 
the commonly used expressions for the lepton asymmetry, which depend on the total CP 
asymmetry and one single efficiency factor, fail to reproduce the correct lepton asymmetry 
in a large number of cases. As an application, we also present in this paper a detailed 
analysis of flavour effects on lepton asymmetries for a two right-handed neutrino model. 
Explicit examples in which sizeable enhancements can be obtained are also given. 

The paper is organized as follows. In Section 2, we review the conventional flavoured- 
blind computation of the baryon asymmetry, and present some analytic approximations 
that will be used later. In Section 3, we introduce the Boltzmann equations we will solve, 
which differ by the inclusion of flavour and of CP violation in AL = 1 processes. The 
following section provides a list of rules and expressions to apply in order to obtain an 
estimate of the baryon asymmetry which includes flavour effects. Section 5 contains the 
analysis in the context of two right-handed neutrino model which make manifest the dif- 
ference between the results when flavours are and are not included. In section 6 different 
textures for the neutrino Yukawa coupling matrix and their implications are explored. In 
Section 7 we discuss the special case in which there is no CP violation in the right-handed 
neutrino sector, and finally in Section 8 we draw our conclusions. 



2 The conventional computation of the baryon asym- 
metry 

This section introduces notation and reviews the calculation of the lepton asymmetry when 
the charged lepton Yukawa couplings are neglected. As we shall see, the commonly used 
formulae for the final lepton asymmetry, which we report in this section, may not be 
appropriate once flavours are considered. 

Our starting point is the Lagrangian of the Standard Model (SM) with the addition of 
three right-handed neutrinos Ni [i = 1, 2, 3) with heavy Majorana masses M3 > M2 > Mi 
and Yukawa couplings Xia- Working in the basis in which the Yukawa couplings for the 
charged leptons are diagonal, the Lagrangian reads 

C = CsM+ (^^N^ + KaN^i^ H + h^H'' enJo. + h.c.^ . (1) 

Here and e^a indicate the lepton doublet and singlet with flavour [a = e, /i, r) respec- 
tively, and H is the Higgs doublet whose neutral component has a vacuum expectation 
value equal to u = 175 GeV. 

We assume that right-handed neutrinos are hierarchical, M2,3 ^ Mi so that studying 
the evolution of the number density of A'^i suffices. The final amount of {B — C) asymmetry 
can be parametrized as Yj^^c = f^B-c/s, where s = 2Tc'^g^T^ /45 is the entropy density and 
counts the effective number of spin-degrees of freedom in thermal equilibrium [g^, = 217/2 in 
the SM with a single generation of right-handed neutrinos). After reprocessing by sphaleron 
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transitions, the baryon asymmetry is related to the L asymmetry by [15] 

where nji^ is the number of Higgs doublets, and uq the number of fermion generations (in 
equilibrium) . 

One defines the CP asymmetry generated by N\ decays as 



UO W R-rr 2-^ 



Im [iXX%] (M, 



[AAt; 
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where the wavefunction plus vertex contributions are included in [T6 



1 — X 



+ 1 - (1 + x)ln 



X 



X 



2:>1 



(4) 



Notice, in particular, that ei denotes the CP asymmetry in the total number (the trace) of 
flavours. 

Besides the CP parameter ei, the final baryon asymmetry depends on a single wash-out 
parameter. 



K = 



mi 
rh* 



(5) 



HiMi) 

where H{Mi) denotes the value of the Hubble rate evaluated at a temperature T = Mi 
{m* ~ 10"^ eV) and 



nil 



Ml 



(6) 



is proportional to the total decay rate of the right-handed neutrino A^i. 

By defining the variable z = Mi/T, the Boltzmann equations for the lepton asymmetry 
Yc, and the right-handed neutrino number density Yn^^ (both normalised to the entropy s), 
may be written in a compact form as 



- Y^l 



dz 



dYr 



dz sH{Mi 



sH{Mi) 
Yn 



ilD + 7al=i) 



Yj 



1 - 



dY^r^ 



dz 



k - I] ei7D - 



Yr 



Y, 



EQ 



ilD + lAL=l + 7AL=2) 



(7) 



(8) 



The processes taken into account in these equations are decays and inverse decays with 
rate 7d, AL = 1 scatterings such as {qf^ Ni), and AL = 2 processes mediated by heavy 
neutrinos. The first three modify the abundance of the lightest right-handed neutrinos. 
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The AL = 2 scatterings mediated by A^2,3 are neglected in our analysis for simplicit}0. The 
various 7 are thermally averaged rates, including all contributions summed over flavour 
(s, t channel interference etc); explicit expressions can be found in the literature (see for 
example [SIE]). Notice that in this "usual" analysis, AL = 1 scattering contributes to the 
creation of A'^i's and not to the production of a lepton asymmetry, only to the washout 
0. This is a minor point in the single flavour analysis; it is more relevant when flavour is 
included, and will be discussed in the following two sections. 

Approximate analytic solutions for Yc and A^r-^ = Y^v^ — Y^^, which reproduce the 
numerical plots [51 E], can be obtained from simplified equations [HIE]- Calculating in zero 
temperature field theory for simplicit}!^, one obtains 



The Boltzmann equations can be approximated 

A^, = -zK^f\{z)A^,-Y^^', (10) 

Y'^ = e,Kz^^A^,--z'KK,{z)Uz)Yc (11) 

where Ki and K2 are modified Bessel functions of the second kind. The function fi{z) 
accounts for the presence of AL = 1 scatterings [3, [6], and f2{z) accounts for scatterings in 
the washout term of the asymmetry. They can be approximated [U], in interesting limits, 

as 

, , f 1 for ^ > 1 

/iW-i N^M for.<l ^^^^ 



and 



1 for ^ > 1 

87r2«222 



Uz)^{^_^ for.<l, 



where ^2^2 = Ks/K ~ 0.1 parametrizes the strength of the AL = 1 scatterings and 
ax = 4/3(2) for the weak (strong) wash out case. A good approximation to the rate 
K z{Ki{z) / K2{z)) fi{z) is given by the function [Kg + Kz) [5l [6] while the wash out term 
— {l/4:)z^KKi{z)f2{z)Yc is well approximated at small z by —aKKgYc- 

In the strong wash-out regime, the parameter K ^ 1 and the right handed neutrinos 
iVi's are nearly in thermal equilibrium. Under these circumstances, one can set A'^y^ ~ 
and A7V1 — {zK2/4:g^K). Exploiting a saddle-point approximation in Eq. ffTTl) one easily 
reproduces the fit to the numerical results [3, E] 

ei /0.55 X 10^3 eV\^-^^ , 
Yc ~ 0.3 - . . (14) 



9* \ rrii 



^Sec,e.g. the Appendix of [13]. We discuss later the restrictions this imphes. 

^We thank A Strumia, A. Pilaftsis, G. Giudice and E. Nardi for useful conversations about this point. 
^Significant finite temperature corrections were found in [S], which have 0(1) effects on the final asym- 
metry. 
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In the opposite weak wash-out regime, assuming that no right-handed neutrinos are 
initially present in the plasma, there could be a cancellation in the final lepton asymmetry 
between the (anti-) asymmetry generated in Ni production, and the lepton asymmetry 
produced as the A''i decay. However, this cancellation does not occur, in Eqs. and 
(IE]), because CP violation in the AL = 1 scatterings is not included. AL = 1 processes 
contribute significantly to the production of right-handed neutrinos, without making any 
associated (anti-) lepton asymmetry, and the A'^i later produce a lepton asymmetry in decay. 
The number of A^^i produced is oc K, and the final lepton asymmetry can be approximated 

m 

Yc=^0.3'-l( „„ "L. \. (15) 



V3.3 X 10-3 eV 

Notice that the final baryon asymmetry in the 'one fiavour approximation' depends always 
upon the trace of the CP asymmetries over fiavours, ei, times a function of the trace over 
fiavours of the decay rate of the right-handed neutrinos, K. This is due to the fact that 
the inverse decay term in Eq. is proportional to the trace over the fiavours of the 
lepton asymmetry times the trace over the fiavours of the decay rate of A'^i's. The reader 
is invited to remember this point in the following when we explain why the 'one fiavour 
approximation' fails to predict the exact baryon asymmetry. 



3 Including flavours and CP violation in scattering 

In this section, we introduce the Boltzmann equations for individual fiavour asymmetries, 
[8]. We define Y^a to be the lepton asymmetry in fiavour a, where the a are the lepton mass 
eigenstates at the temperature of leptogenesis. As discussed in [13], the Yaa are the diagonal 
elements of a matrix [Y] in fiavour space, whose trace is the total lepton asymmetry. In 
this paper, the off-diagonal elements are neglected 0. 

The mass eigenstates for the particles in the Boltzmann equations(BE) are determined 
by the interactions which are fast compared to those processes included in the BE. The 
interaction rate for Yukawa coupling ha can be estimated as [17] 

Ta^5xlO~'^hlT, (16) 

so interactions involving the tau (mu) Yukawa coupling are out-of-equilibrium in the 
primeval plasma if T > lO^^ GeV (T > 10^ GeV) |. Thermal leptogenesis takes place at 
temperatures on the order of Mi, and the asymmetry is generated when the rates ^ H, 
so we conclude that the r (/i) lepton doublet is a distinguishable mass eigenstate, for the 
purposes of leptogenesis, at T < 10^^(10^) GeV. 

^See [13] for a discussion. The equations of motion for the matrix [Y] are more comphcated than the 
Boltzmann equations, but at most temperatures are equivalent to Boltzmann equations written in the mass 
eigenstate basis of the leptons in the plasma. The off-diagonal elements of [Y] could have some effect on 
the lepton asymmetry, if leptogenesis takes place just as a charged lepton Yukawa coupling is coming into 
equilibrium (so the mass eigenstate basis is changing). 

^The electron Yukawa coupling mediates interactions relevant in the early Universe only for temperatures 
beneath ^ 10^ GeV and can be safely disregarded. 
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The Boltzmann equations that we will use in this paper, for the flavour asymmetries 
Yaa, are listed below. They differ from those of [13] in two respects. 

First, we have neglected the off-diagonal terms of the matrix [Y]. The second and 
most significant difference is that we have included CP violation in the AL = 1 scattering 
rate, which will give Yl cx for weak washout (instead of K in eq. ( |T5|) ). That is, the 
function 7al=i, which appears in the A'^i creation term, also now appears in the first term 
of equation (fT71) . which describes the production of the lepton flavour asymmetry. Later in 
this section we will calculate the CP asymmetry in scattering, and show that it gives the 
same e as decay and inverse decay, in the limit of hierarchical right-handed neutrinos. As 
in [13] , we continue to neglect the non- resonant contribution to AL = 2 scatterings, and its 
associated flavour effects [Sj. At the end of section 4.1.2, we discuss the parameter range 
where this is acceptable. 

Equation (JTj) for the A'^i number density remains unchanged, and the equation for the 
flavoured lepton asymmetry is 



dY^ 



- 1 e (7z) + 7AL=l) - --Ro (7Z3 + 7al=i; 



^EQ ^ p ' '^^=^1 ^EQ 



17) 



dz sH{Mi] 

where there is no sum over a in the last term of equation (1171) (or ffTH]) ). 

The rates [H] and asymmetries are calculated in zero temperature fleld theory, and 
include processes mediated by the neutrino and top Yukawa couplings. This is simple, 
and parametrically consistent. However, flnite temperature and gauge corrections can be 
signiflcant [5], and should in principle be included. 

To obtain analytic solutions we simplify this, with the approximations introduced in 
the previous section, to 



Y^a = eaaKz^^hiz)AM,-]z'Ki{z)f2{z)KaaYaa (1^ 

A 2 (2) 4 



AW = -^K^^Mz)A^,-Y^^', (19) 

where 

Kaa parametrizes the decay rate of A^i to the a-th flavour, and the trace K^a, coincides 
with the K parameter deflned in the previous section, see Eq. ([5]). 

Notice in particular that the dynamics of the right-handed neutrinos is always set by 
the total K. 

The CP asymmetry in the a-th flavour is and is normalised by the total decay rate 

= (21) 

' -Im LJ^xJ (22) 



(167r)|AAt]ii [ '» V 
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where the second hne is in the hmit of hierarchical Nj, and m = U*DmU'^ = v'^X^ M~^X is 
the hght neutrino mass matrix. If rrimax is the heaviest hght neutrino mass (= matm for the 
non-degenerate case) and we define e"^^^ = 3Aml^^Mi/{8nv'^mmax) [IH], then the flavour 
dependent CP asymmetries are bounded by 



3Mimmax I _ max ^max h^aa /r,QN 

w V"r-' '^^^^ ^^^^ 

so the maximum CP asymmetry in a given flavour is unsuppressed for degenerate light 
neutrinos [13], but decreases as the square root of the branching ratio to that flavour 

The CP asymmetry eaa can be written in terms of the diagonal matrix of the light 
neutrino mass eigenvalues m = Diag(mi, m2, ma), the diagonal matrix of the the right 
handed neutrino masses M = Diag(Mi, M2, M3) and an orthogonal complex matrix R = 
w M~^/^ A f/m~^/^ [in], where U is the leptonic mixing matrix, which ensures that the 
correct low-energy parameters are obtained 



3Mi Im (E,3p"^y^"^p^^^a/3f^"P^/3i^pi 



As noted in [T3], for a real R matrix, the individual CP asymmetries niay not vanish 
because of the presence of CP violation in the U matrix. On the contrary, the total CP 
asymmetry ei = J2a vanishes. 



3.1 The CP asymmetry in AL = 1 scattering 

We now wish to show that the CP asymmetry in scattering processes such as [qf^ — > Nl^) 
is the same as in decays and inverse decays!^. This result was found in [HUSO], for the case of 
resonant leptogenesis. CP violation in scattering is usually neglected in thermal leptogenesis 
[SIE], because observed neutrino masses favour the strong washout regime K > 1. In strong 
washout, any contribution to the lepton asymmetry by scattering processes during A'^i 
production is rapidly washed out, so the CP violation in these processes can be neglected, 
for lack of observable consequences. However, we wish to include this CP violation, because 
washout in one flavour could be small, even though K ^ 1 and we wish to correctly include 
the contribution of weakly washed out lepton flavours to the final lepton asymmetry. 

For simplicity, we work at zero temperature, in the limit of hierachical right-handed 
neutrinos. This means we calculate in an effective field theory with particle content of the 
SM -t-A^^i, and the effects of the heavier N2 and N^, appear in a dimension five operator 
{HLa){HLp). For computing one-loop CP violating effects involving A^i, we can take the 
coefficient of this operator oc [mi,]a/3 / . 

We define the CP asymmetries in AL = 1 scattering (mediated by s and t-channel 
^We thank E. Nardi for discussions of his work in progress. 
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Higgs boson exchange) as 



a + o 

(T{qN -^FL^)-cx{qN -^t'L^) a{qL^ FN) - cx{qL^ ^ t'N) 

e, = = (2b) 

* a + a a + a ^ ' 

where barred fields are the antiparticles. The initial state density factors cancel in the ratio, 
so the cross-sections a, a can be replaced by the matrix elements squared integrated 
over final state phase space J (ill. If the tree + loop matrix element is separated into a 
coupling constant part c and an amplitude A: 

M = c'A' + c^A^ , (27) 

where the matrix element for the CP conjugate process is Ai = c^*A^ + c^*A^, then the CP 
asymmetry can be written 

_ 21m{c*c'*} / Im{^U'*}rfn 

The loop amplitude has an imaginary part when there are branch cuts corresponding to 
intermediate on-shell particles, which can arise here in a bubble on the TV line at the NHLa 
vertex, e.g. for s-channel Higgs exchange: 

lm{A'A^*} = A'if'q NL^) j A'*{t^q L^L'f^H')dU'A'*{L'pH' N) . (29) 

H' and L'^ are the (assumed massless) intermediate on-shell particles, and dU' is the inte- 
gration over their phase space. 

In the scattering process c* = htXl^ and = 3ht\ii3['m*]pa/ v"^ , where ht is the top 
Yukawa coupling, so the complex coupling constant combination in iaai^H — > NLa) is 
clearly the same as in eaa of eq. ( l22ll^ . To obtain the amplitude ratio (the second fraction 
in eqn (128|) ). we take, for instance A^{N HL'^) = uePi,U]\f, and after straightforward 
spin sums, one sees that it is the same for scattering and decay, so e"" = = e°"^. 

The equality (1261) . and the result that the s and t-channel scattering asymmetries are 
equal to the decay asymmetry, may be an artifact of our effective field theory calculation, 
where there is no momentum exchanged on the iV2,3 line in the loop. 

In the finite temperature calculation of [5] , the Higgs boson has a large thermal mass due 
to its interactions with the top quarks. At T ^ Mi it can therefore decay to NL, and ref. 
[5] found a CP asymmetry in this decay. Some part of the N production that is included in 
zero temperature scattering computations, is resummed in the finite temperature H decays, 
so it is consistent that in both processes, a lepton asymmetry can be generated. 



^°In both cases, there is an overall factor of three, from the weak SU{2) index contractions. This can be 
seen by reinstating the A^2,3 propagators: the charged and neutral component of the intermediate H' and 
L' contribute in the iVi wave-function correction, giving a 2, but only the charged or the neutral H' and 
L' appear in the vertex correction. 
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Checking CPT and unitarity 



We would like to verify that the CP violating matrix elements M. for AL = 1 scattering 
processes, satisfy the constraints following from CPT invariance and S'-matrix unitarity. 
See e.g. [21] for useful related discussions. For i any state, and i its CP conjugate, CPT 
and unitarity imply 



E \MiJ = \MCj ^^)\' = Y1 ^ ^Or, (30) 

j j j 

where {j} is the set of accessible states. We want to check that this is consistent with 
having a CP asymmetry in scattering processes like qf^ Ni. 

First, as a warm-up, let us consider the case of inverse decays Hi —>■ N in the single 
flavour approximation. It is well-known that there is a CP asymmetry here, and this will 
show how to apply CPT and unitarity bounds with unstable particles (N) in the final state. 

Suppose that 

I 12 I A/l |2 

mm - iv)p = ^(1 + 6) , \Mim^N)\' ^ ^(1 - 6) , (31) 

where |A1oP is order and |A^ope is order A^. 

At order A^, the CPT + unitarity constraint is verified, but at order A^, \A4{Hi -^■ 
A^)P 7^ \M{Hi — > A^)P, because there can be additional final states: HI and H£. As 
in the case of the Boltzman Equations, when we include these 2 2 processes, we must 
subtract out the real intermediate A'^i, because we have already included this by treating 
the Ni as a final state particle. So at order A^: 

\M{Hi^X)\^ = \M{Hi-^ N)\^ + \M{Hi^ m)\^ -\M^^^{H£^ Hi)\^ 
+ \M{Hi-^Hi)\^-\M^^^iHe^He)\^ 

"■^°"'a + e) + |Al(i/£^M)r-^(l + e)il±il 



2 

+ \M{Hi^Hi)\ 



2 \Mo\\, , Al-e) 



2 ' ' 2 

= \M{Hi^Hi)\^ + \M{Hi^Hi)\^ + --- (32) 

where X is all possible final states, and the on-shell intermediate A'^i has been subtracted 
in the narrow width approximation [S]. Repeating for {Hi — X) will give the same 
rate, because at O(A^), \M{Hi Hi)\^ = \M{Hi Hi)\^, and \MiH£ H£)\^ = 
\Ai{Hi Hi)\'^. So CPT and unitarity are satisfied for inverse decays, as they ought to 
be. 

CPT and unitarity are realised in the scattering process {qf^ —>■ Ni^), in a similar way 
to inverse decays. CPT and unitarity should hold order by order in perturbation theory, 
so we work at order X^X^h^, and define 

\M{qf^NQ\^ = |A^,|2(l + e"-), (33) 
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where \Ms\'^ oc A^/i^, and |7Wspe oc X'^X'^hf. At order X'^X'^hf, we should also include 
various tree diagrams without in the final state. Following the inverse decay discussion, 
one can write 

+ \M{qt^^lHi^)\'-\Ms\\l + en^^-^ 

= J2 0-^« + ^ , (34) 

where, in the narrow width approximation, 

IM^^^iqt" IpHQ\^ = \M{qt^ Ni^)\^ x BR{N Hip) (35) 

In eq. flH^ . the CP asymmetry e"" has disappeared, so if we repeat the calculation for 
the CP conjugate initial state qf^, we should obtain the same result, verifying that a CP 
asymmetry in qf^ — s> N^a is consistent with CPT and unitarity. Furthermore, eq. (!34|) is 
reassuring, because the unstable state has disappeared. There is no CP violation in the 
total rate for qf^ asymptotic (stable) final states, but CP violation in the partial rate to 
the unstable A^ is possible. This can be relevant to the final value of the baryon asymmetry 
when some of the lepton fiavours are weakly washed out (see section l4.1.2p . 



4 Approximate formulae for the baryon asymmetry 

In this section we present analytical formulae for the final baryon asymmetry in the case 
in which fiavours are taken into account. We can have different possible cases according to 
which interaction mediated by the charged Yukawa couplings is in equilibrium, [HI [H] . 

4.1 ji and r Yukawa couplings in equilibrium: Mi < 10^ GeV 

The /i and r doublet leptons (and by default, the electron doublet) will be mass eigenstates 
at the temperature of leptogenesis when the mass of the lightest right-handed neutrino Mi 
is smaller than about 10^ GeV. The Boltzmann equation for the diagonal entry Y^a reads 
(no summation over the index a) 

YL = ^caKz ^l^fi{z)AM, - ]z'Ki{z) f2{z) K^aYaa • (36) 

K2{z) 4 

Let us now solve analytically Eq. (l36l) according to the magnitude of the various Kaa- 
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4.1.1 Strong wash-out regime for all flavours 



In such a case all the K^a ^ 1- The right-handed neutrinos A'^i's are nearly in thermal 
equilibrium. Under these circumstances, one can set A^^ ~ and A^Vj — {ZK2/ '^g^.K). 
The lepton asymmetry for the flavour a is given by 

Jo 

Using the steepest descent method to evaluate the integral, one flnds that it gets the major 
contribution at 'z such that 2 = log Kaa + (5 In 2/2) when inverse decays become ineflicient. 
The lepton asymmetry in the flavour a becomes 

e„„ /0.55 X 10-3 eV^ ^'^^ 



y,,^o.3— r-"'r" . (38) 

To get convinced that this result differs from the one usually considered in the literature, 
let us take the total lepton asymmetry = F^a 

y,^Vo.3f^f°-^S'<l°"^Vy'°. (39) 

^ g* \ fllaa J 

It does not coincide with total lepton asymmetry result (fT^ in the strong wash-out regime. 
Indeed, the total lepton asymmetry flH^ is the sum of the eaa, each weighted by the wash- 
out factor Kaa and not the sum of the divided by the sum of the Kaa- Eqs. ffT^ and 
( !39|) coincide only if one family is dominating the contribution to the total CP asymmetry 
and the corresponding wash-out factor is the tiniest. 



4.1.2 Weak wash-out regime for all flavours 

In this case all the Kaa ^ 1- We assume that right-handed neutrinos are not initially 
present in the plasma, but they are generated by inverse decays and scatterings. As ex- 
plained in Sec. 2, the equation of motion for Yati is well approximated by 

= -{K, + Kz) (y^v. - Y^f) , (40) 

We split the solution in two pieces. Let us deflne z-^q the value of z at which Y/Vi(^eq) = 
Y^'^{z-^q). This value has to be found a posteriori. For z <^ 2;eq, we may suppose that 
Yn^ < Yn? and Eq. (00]) is solved by 



Y^^{z)c:, r dz'{Ks + Kz')Y^^ = r dz'{Ks + Kz'){zTK2{z' 
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With Ii and I2 integral involving the modified Bessel functions : 

hiz)= [ x^K2ix)dxc^ fiz) + z^K2iz) (42) 
Jo 



(44) 



where [6] 

f{z) = — 7r,c = 0.7 (43) 

((97r)- + (2z3)c)Vc' 

The integral I2 is well known, and equals 

/2(z) = / x^K2{x)dx = 8- z^K3{z) 
Jo 

Therefore 

YnM ^ ^ + ^'^2(^)) + 8 - z'K,{z)^ (45) 

As expected for weak washout, we find that the maximum number density of A^^i is 
proportional to K (recall Ks oc K). 

Let us now compute the value of zeq- We expect it to be ^ 1 and we therefore 
approximate, up to C(z~^/^) : K2{z) ^ Ks{z) ~ ^^z~^/'^e~^. Imposing Y^_^{zeq) = 

^Wi^(^eq), we find 

ZEQ ^ ^ In ZEQ - In (^-^^K + ^^K,^ . (46) 

This solution is a good approximation to the real value for i^T <^ 1. 
For z > zeq, we have {Kg + Kz) ~ Kz and 

r^,(^)^r^Q(^EQ)e^/^(^-Q-^'). (47) 

We have included CP violation in AL = 1 scattering, unlike the usual analysis, so we 
expect our solution for Y^" to have a different scaling with K than eqn ffT^ . The reason 
is as discussed in [SI E] : if (^P in scattering is neglected, then Yjy oc K, and the A^i decay 
out of equilibrium, so one expects Y£'" oc Ke°"^. However, if (fP in Ni production (~ 
scattering) is included, and washout is neglected, then the equations for Y^^ and Y^" are 
identical, so Y^'^{z — * 00) vanishes. That is, for every |l/e""| A''i's that are created, be 
it by inverse decay or scattering, an (anti)-lepton a is produced. This (anti-) asymmetry 
will approximately cancel against the lepton asymmetry generated later on, when the A^^i 
decay. However the cancellation will be imperfect, because the anti- asymmetry has more 
time to be washed out, so the final asymmetry should scale as KKaa- After integrating by 
parts, this is what we find for the asymmetry in the fiavour a, which is given by 

Y^^ - e^a rfz'F;vi(^')^?aa(^')e~^^°°'''''"'^^''^ ^7aa(^) = 7;2'i^aai^i(^)/2(^), 

Jo 4 

^ 1.5^f ^ ]( ^ V (48) 

V3-3 X 10-3 eV; V3-3 X 10-3 eV; ^ ' 
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We have checked numerically that this analytical formula fits the numerical results to a 
30%. 

Our findings hold provided that the non-resonant AL = 2 scattering rates, in particular 
those mediated by the N2 and heavy neutrinos, are slower than decays and AL = 1 
scatterings when most of the asymmetry is generated. We estimate that this applies when 

This means that K^a should be larger than 10~^. 



4.1.3 Strong wash-out for some flavours and either weak or mild wash-out for 
others 

In the case in which Kaa ^ 1 for some flavour a, but Kpp -C 1 for some other flavour /3, it 
is impossible to match this case with any of the cases discussed in the section for the 'one- 
flavour' approximation. Because some flavours a strongly interact with the right-handed 
neutrinos, we may set K = Kaa ^ 1- Thus, right-handed neutrinos are brought 
to thermal equilibrium by inverse decays and by AL = 1 scatterings to an abundance 
approximately given by 

Y^,{z):^Y^^{l-e'^^^-^^y') . (50) 
The lepton asymmetry in flavour /3 is then given by 



^ Oa'-^ ) . (51) 

\3.3 X 10-3 eV^ ^ ' 

Again, we have checked that this formula fits the numerical results to about 30 %. It has the 
dependence on K^/j expected from our inclusion of (J^P in scattering. The anti-asymmetry 
~ —epps/g^ created during A''i production has the leisure to be partially destroyed by 
processes violating (3 lepton number. The amount of this reduction is controlled by K^j^, 
so one expects the final [3 asymmetry oc ej3pKpp, as in eq. ( jSTl) . Again, this result applies 

for Kaa > 10-1 

In the case in which Kaa ^ 1 for some flavour a and the flavour (3 suffers a mild 
{Kjjp ~ 1) wash-out, the final asymmetry in the flavour (3 is fitted within 30 % by the 
formula 

'^^-^I^U.25xlO-3evJ +1 ) ) ■ ^^'^ 

This formula reduces to Eqs. (|38l) and (jSTl) for the weak and strong wash out for the flavour 
P, respectively. 
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4.1.4 Recipe to go from the flavour asymmetries to the baryon asymmetry 

The final baryon asymmetry is given by [110 

12 



37 



E , (53) 



where the asymmetries Y/\^, (Aq = {B/3 — La)) are conserved by the sphaleron and other 
Standard Model interactions. To provide an analytical formulae for the baryon asymmetry, 
in the case in which Mi ^ 10^ GeV, it suffices to obtain the asymmetries from our 
formulae for the lepton ffavour asymmetries Yaa- 

In this paper, we approximate = Yaai'^^aa ^aa^aa)- That is, in our approximate 
solutions for Yaa-, we replace fhaa by Aaa'fnaa (no sum on a), where the matrix A [8] is 
given in eqns (15^ and ( l56l) . This approximation should be adequate for transforming the 
ffavour asymmetries, estimated ffavour by ffavour, into a baryon asymmetry. To see this, 
we review the discussion of [H] : 

The interactions which are much faster than H , such as gauge interactions, sphalerons, 
and some Yukawas, can be approximately taken into account by imposing chemical equi- 
librium conditions(see [211 E] ^ more careful analysis). This allows to express the 
asymmetries in all particle species as linear combinations of the asymmetries in conserved 
quantum numbers. These interactions-in-equilibrium respect the three ffavoured asymme- 
tries = B /3 — La, which are changed by the interactions of N. So it is the Boltzmann 
Equation for Y^^ that is relevant for leptogenesis, and it is clearly [S] the right-hand-side 
of eqn 



dY, 



Ac 



dz sH{Mi] 



1 e (7b + 7al=i) - — EQ {Id + 7al=i) 



Yaa on the right-hand-side can be re-expressed as Yaa = X]/3^a/3^A^- An ia produced or 
destroyed in the plasma changes 5/3 — L^, but not all the B/3 — La is in Yaa, so washout 
is reduced and the ffavoured Boltzmann Equations are coupled. 

In the temperature range where the /i Yukawa is in equilibrium, the A matrix is given 
as@ 

-151/179 20/179 20/179 \ f Y^\ 
25/358 -344/537 14/537 Y^^ . (54) 

25/358 14/537 -344/537 / \ Ya, J 

The off-diagonal elements of A couple the ffavoured Boltzmann Equations. We neglect 
this effect, estimating that it has little effect on the ffnal baryon asymmetry: since the 
off-diagonals of A are small, AapY^^^ (for f] ^ a, not summed) will only make a signiffcant 
contribution to the washout of when 1^^ ^ ^Aq • In this case, the contribution of Y^^ 



^^if the sphalerons freeze out before the eletroweak phase transition. 12/37 would be replaced by another 
fractiongS] of order 1/3, such as^ 28/79. 

^^This disagrees with our previous version[T3j, where we had taken the u and d Yukawas in equilibrium, 
and differs from [8] in that we include strong sphalerons. Our results agree with [24} I14j. 
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to the baryon asymmetry is small, so it is of little numerical importance that our formula 
mis-estimated the washout in flavour a. 

Thus the recipe to compute the baryon asymmetry in the case in which Mi ^10^ GeV 
is the following: 

1) Compute the K^a parameters for the three flavours to see which of the three different 
previously discussed cases applies; 

2) compute the CP asymmetry in each flavour, and approximate the asymmetry 
using the relevant formualae for Y^a among equations (139|) . ( l48i) . ( l5Ti) . hut with fhaa 

replaced by \Aaa\'fnaa (no sum on a). 

3) compute the baryon asymmetry 



^ 37 



151 . 



344 . 



344 . 



(55) 



Figs. [H [2] and [3] illustrate the differences between the final baryon asymmetry with 
flavours accounted for and the result obtained within the one-flavour approximation for 
some values of the wash-out parameters and CP asymmetries. 



10-'- 




Figure 1: The total baryon asymmetry in the two flavour calculation (upper) and within the 
one-flavour approximation (lower) as a function of z. The chosen parameters are Krr = 10, 
Kf,f, = 30, Kee = 30, err = 2.5 X 10~^ e^^ = -2 x 10~^ eee = lO"'^ and Mi = 10^° GeV. See the 
second Note Added at the beginning of the manuscript. 
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Figure 2: The total baryon asymmetry in tlie two flavour calculation (upper) and within the one- 
flavour approximation (lower) as a function of z. The chosen parameters are K^-t = 4.5 x 10"^, 

= 10-'^ and Ml = 10^° GeV. 
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Note Added at the beginninf) of the manuscript. 
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Figure 3: The total baryon asymmetry in the two flavour calculation (upper) and within the 
one-flavour approximation (lower) as a function of z. The chosen parameters are Ktt = 10, 
Kf,f, = 30, Kee = 10~2, err = 2.5 X 10-^ e^,^ = -2 X 10"^ eee = lO-'^ and Mi = 10^° GeV. See 
the second Note Added at the beginning of the manuscript. 

4.2 Only the tau Yukawa coupling in equilibrium: 

(10^ < Ml < 1012) GeV 

This case is realized when the mass of the hghtest right-handed neutrino Mi is larger than 
about 10^ GeV, but smaller than about 10^^ GeVo. Interactions mediated by hr are in 
equilibrium, but not those mediated by /i^. The off-diagonal entries of the matrix Yrp 
all vanish and the muon and electron asymmetries are indistinguishable. The problem 

^•^In the case in which Mi is around 10^ GeV ofF-diagonal terms may be relevant. However, they are 
quickly damped away as soon as Mi becomes larger than 10^ GeV |13j . 
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of finding tlie total baryon asymmetry reduces to a case of two flavours, the lepton i^, 
and £2, the non-r components of the lepton into which Ni decays. At tree level, £2 = 

Sa=e /X -^la^a/ \ )lia=e p. I"^!"! ) • therefore define two CP asymmetries, err and 

£2 = fee + e^t/i, and the corresponding wash-out parameters K^^ and K2 = K^e + K^j^^j^ for 
the two asymmetries Y^-t and Y2 = Y^e + Y^^. One then solves this set of equations as 
we described in the previous subsection depending upon the magnitude of the wash-out 
parameters K^-r and K2. 

The recipe to compute the baryon asymmetry in the case in which (10^ ^ Mi ^ 10^^) 
GeV is the following: 

1) Compute the K^-t and K2 parameters to see whether they are either both larger than 
unity, or both smaller than unity, or one larger and the other smaller than unity; 

2) compute the CP asymmetry in each of the two flavours, and the 5/3 — asym- 
metry Iac. using the relevant formulae among equations (IHUl) . (HHj) . (1^ . replacinqin these 
formulae rfipp by (no sum on (3), where in this case the A matrixjS] isl 



Y2 \ _ ( -417/589 120/589 \ / Fa^ 
J ~ \ 30/589 -390/589 ) I Fa, 



(56) 



3) compute the baryon asymmetry 



(^2(62, 4177712/589) + i;.(e.., 390m.,/589)) . (57) 



In this case, the first row, second colomn element of the matrix ( l56l) is ~ 30% of the 
diagonal elements, so the approximation of neglecting off-diagonals is less good than in the 
T < 10^ GeV case. Again, our results apply if AL = 2 scatterings are negligible, see eq. 



4.3 No charged Yukawa couplings in equilibrum: Mi ^ 10^^ GeV 

This case is realized when the mass of the lightest right-handed neutrino Mi is larger 
than about 10^^ GeV. Interactions mediated by all charged lepton Yukawa couplings 
are out of equilibrium and all flavours are indistinguishable. The problem of find- 
ing the total baryon asymmetry reduces to a case of one flavour, the lepton £1 = 

Zla=e,M,r -^la^a/ ( I]a=e,AJ,T l"^!"! ) ' ^"^^ ^au therefore define a single CP asymmetry, 
ei = e^Q, and the corresponding wash out parameter K = ^aa- The corresponding 
Boltzmann equations read as those in Eq. fllip . Only in this case the commonly used for- 
mulae reported in Section 2 hold with the eventual inclusion of the effects from the AL = 2 
scatterings. 

^"'This differs slightly from results in [SKU] because we consider the case where the c Yukawa and strong 
sphalerons are in equilibrium 
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5 General analysis of leptogenesis in the two right- 
handed neutrino model 

In this section we will concentrate on the see-saw model with two right-handed neutrinos 
(2RHN). Oscillation experiments indicate that two new mass scales have to be introduced, 
in order to account for the solar and atmospheric mass splittings. These two mass scales 
could be associated to the masses of two right-handed neutrinos, therefore a see-saw model 
with just two right-handed neutrinos can already accommodate all the observations. An 
additional motivation to consider the two right-handed neutrino model is that it corre- 
sponds to some interesting limits of the complete see-saw model with three right-handed 
neutrinos, namely when the mass of the heaviest right-handed neutrino is much larger 
than the masses of the other two, or when the Yukawa couplings for the first generation of 
right-handed neutrinos are much smaller than for the other two generations. 

The two right-handed neutrino model depends on many less parameters than the com- 
plete see-saw model and the theoretical analysis of leptogenesis becomes much more man- 
ageable, while preserving the key features of the model with three right-handed neutrinos. 
In the basis where the charged lepton Yukawa coupling and the right-handed mass ma- 
trices are diagonal, the model is defined at high energies by a 2 x 3 Yukawa matrix and 
two right-handed neutrino masses. Mi and M2. This amounts to eight moduli and three 
phases. On the other hand, at low energies the neutrino mass matrix is defined by five 
moduli (two masses and three mixing angles) plus two phases (the Dirac phase and the 
Majorana phase). In particular, since the mass matrix is rank 2 in this model, the lightest 
neutrino mass eigenvalue automatically vanishes and only two possible spectra may arise: 

• Normal hierarchy: mi = , m2 = y^Am^ , = y^Amf^ 

• Inverted hierarchy: rris = 0, rrii = a/ Am^^„ — Aml^i , m2 = a/Attt^^ 

The only Majorana phase corresponds to the phase difference between the two non- 
vanishing mass eigenvalues. Therefore, the number of unmeasurable parameters in the 
2RHN model is reduced to three moduli and one phase. 

The most general Yukawa coupling compatible with the low energy data is given by: 



where m = Diag(mi, m2, ^3) is the diagonal matrix of the light neutrino mass eigenvalues 
(which has mi = for the normal hierarchy and rris = for the inverted hierarchy), 
M = Diag(Mi,M2) is the diagonal matrix of the right handed neutrino masses, U is the 
leptonic mixing matrix, and R is an orthogonal matrix that in the 2RHN model has the 
following structure [26] 




X = M'/^Rm'/^Uy' 



V 






(normal hierarchy). 



(59) 




(60) 
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with 9 a complex parameter and ^ = ±1 a discrete parameter that accounts for a discrete 
indeterminacy in R. In consequence, the elements of the neutrino Yukawa matrix read: 



Xia = vMi(^/m^cos6' t/*2 + ^ym^sin6' t/*3)/w, 

X2a = ^/M'2i-^2smeu:^ + ^^,coseu:^)/v, (61) 



for the case with normal hierarchy and 



X2a = ^/m'2{- ./T^i sine U*, + ^^/^ cose U*^)/v, (62) 



for the case with inverted hierarchy. The three moduli and the phase that are not de- 
termined by low energy experiments are identified in this parametrization with the two 
right-handed masses Mi and M2, and the complex parameter e. 

Notice that we have included all the low energy phases in the definition of the matrix 
U, i.e. we have written the leptonic mixing matrix in the form U = V Diag(l, e"**^/^, 1), 
where is the Majorana phase and V has the form of the CKM matrix: 

C13C12 C13S12 sise"*"^" 

V = I -C23S12 - S23Sl3Ci2e*'^ C23C12 - S23Si3Si2e*^ S23C13 | , (63) 
S23S12 - £235130126*"^ -S23C12 - C23Si3Si2e*'' C23C13 

SO that the neutrino mass matrix is = ?7*Diag(mi, m2, m3)f/''". It is straightfor- 
ward to check that the Yukawa coupling eq. fl58p indeed satisfies the see-saw formula 
M = X^Biag{Mi\M^^)Xv^. 

The flavour CP asymmetries can be readily computed in terms of low energy data and 
the unmeasurable parameters Mi, M2 and e substituting the expression for the Yukawa 
coupling in eq. (l2T]) . In the limit Mi <C M2 we obtain for the case with normal hierarchy 
the following result: 

eaa - -3-^ ^ T77~~^^ , 7^[("^3l^a3|^ -"^2|f42n Imsin^^ + 

onv 7712! cos^ y| + msl sin e\ 

+C v^m2m3 (777,3 + Re[/*2t/a3 Im sin e cos e + 

+^V"^2"^3("^3 - ^,2) lmU*2Ua3 Re siu e cos e], (64) 

and analogously for the case with inverted hierarchy, with the changes in the labels 3^2 
and 2^1. 

We will analyze numerically the predictions for the baryon asymmetry taking flavour 
properly into account, as described in Section 3, and also for comparison following the 
conventional computation ignoring flavour, as described in Section 2. We will perform this 
analysis both for the case of normal hierarchy and inverted hierarchy, fixing the atmo- 
spheric and solar mass splittings and mixing angles to the values suggested by oscillation 
experiments, Aml^^ ~ 2.2 x 10"^ eV^, Am^„; ^ 8.1 x 10"^ eV^, 6^23 ~ 7r/4 and 6^12 ~ 7r/6, 
respectively [25]. The remaining parameters in the leptonic mixing matrix are fixed to 
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^13 = 0.1, 5 = 7r/4, </) = 7r/3. In general, the results are not very sensitive to the value of 
6*13. On the other hand, the aspect of the plots does depend on the precise value of the 
phases 6 and 0, although our main conclusions remain valid. 

We will show our results for different values of Mi, to cover the possibility that only 
the tau Yukawa interaction is in equilibrium (Mi ^ 10^ GeV) or that the tau and the muon 
Yukawa interactions are in equilibrium (Mi ^ 10^ GeV), and for different values of the 
complex parameter 9, restricting ourselves to the region |^| < 1. We will also fix ,^ = 1, 
although the results for the case ^ = — 1 can be read from our results by changing 6 —6, 
as can be checked from eq. fl64p . 

In FigJH we show the result of the calculation of the baryon asymmetry in the complex 
plane of 6 when the mass spectrum presents a normal hierarchy, following the calcula- 
tion that takes flavour properly into account (left plots) and following the conventional 
calculation ignoring flavour (right plots). In the upper plots we show the results when 
Ml = 10* GeV, so that the tau and muon Yukawa interactions are in equilibrium, whereas 
in the lower plots we take Mi = 10^° GeV, so that only the tau Yukawa interaction is in 
equilibrium. It is apparent from these plots that the proper treatment of flavour in the 
Boltzmann equations is necessary in order to calculate correctly the baryon asymmetry!^ 

The differences between the correct analysis of leptogenesis, taking flavour into account, 
and the conventional analysis are more acute along the axes Im^ = and Re^ = 0, and 
around the values of 6 that correspond to texture zeros in the Yukawa coupling. The 
corresponding values of 6 are indicated in the plots with a red star, when |Ai3| ~ 0, and 
with a blue cross, when |Ai2| — (a texture zero in the (1,1) position appears at |^| > 1). 
The point in the 6 complex plane where \Xia\ — can be derived from eq. flUT]) . being the 
result: 

tan^r - -^\—Jfr- (65) 

V ''''3 '-^aS 

The reason why in the conventional analysis the baryon asymmetry vanishes along the 
axes can be easily understood from the expression of the total CP asymmetry in the 2RHN 
model. For instance, for the case with normal hierarchy: 

3Mi (mo — m^) Imsin^ ^ , , 

ei ~ ^ ^- (66) 

OTTf ^ 7722 I C0S2 6\ + m^l slu^ 6 \ 

Hence, the total CP asymmetry vanishes when Im6' = {i.e. R real) and when Re^ = 
(since in this case sin^ 6 = — sinh^ \6\, and R is again real). This is not necessarily the case 
when flavour is properly taken into account, as can be realized from the expressions of the 
flavour CP asymmetries, eq. fl22l) . 

On the other hand, when the Yukawa coupling presents an approximate texture zero, 
the difference stems mainly from the washing-out of the asymmetry. When the spectrum 



^^In particular, the prediction for the baryon asymmetry in the conventional computation is symmetric 
under Re0 —Kc6 and luiO — > —hn9, and consequently independent of the discrete parameter ^. On the 
other hand, when flavour is taken into account, the parameter ^ indeed plays a role in the computation of 
the baryon asymmetry. 
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Figure 4: Baryon asymmetry in the 2RHN model when the neutrinos present a spectrum with 
normal hierarchy. In the left plots we show the result of the calculation that takes flavour into 
account, whereas in the right plots, the result of the conventional calculation that ignores flavour. 
On the top plots we take Mi = lO^GeV, so that the tau and muon Yukawa couplings are in 
equilibrium, while in the lower plots, Mi = lO^'^GeV, so that only the tau Yukawa coupling is in 
equilibrium. In these plots we have fixed ^13 = 0.1, 5 = 7r/4, cj) = ir/S and the remaining neutrino 
parameters to their favoured experimental values. 

has a normal hierarchy, the CP asymmetry in the a-th flavour is comparable to the total 
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CP asymmetry. However, the smallness of the interaction of the hghtest right-handed 
neutrino with translates into a weak wash-out of the asymmetry, in stark contrast with 
the result of the same analysis following the conventional computation, where the total 
CP asymmetry is necessarily strongly washed-out (recall that in the 2RHN model with 
normal hierarchy rhi > a/ Am^^^, so > 1). As a consequence, when Aia ~ 0, the 
actual prediction for the baryon asymmetry can be around one order of magnitude larger 
than previously believed. Interestingly enough, in realistic models strict texture zeros 
rarely appear; there are usually subleading effects that produce small entries instead of 
strict texture zeros. To be precise, texture zeros normally appear in a basis where the 
charged lepton Yukawa coupling and/or the right-handed neutrino mass matrix are slightly 
non diagonal. Therefore, the diagonalization of these matricial couplings to express the 
Lagrangian as in eq.([T]) will lift the texture zero in the neutrino Yukawa matrix, yielding 
a small entry instead. Finally, even if the texture zero was strict at the high energy scale, 
radiative corrections could lift them. 

The differences between the computation of the baryon asymmetry taking into account 
flavour or not are even more acute in the case of a spectrum with inverted hierarchy, as 
can be realized from Fig. [5l As in the case of the spectrum with normal hierarchy, the 
maximum differences arise along the axes Im6' = and Re^ = 0, and around the values of 
6 that correspond to texture zeros in the first row of the Yukawa coupling, also indicated 
in these plots with a red star and a blue cross 0. The precise value of 9 where |Aiq,| ^ is: 

t&ney - -^\ —jjr- (67) 

Around this value for 6', the difference between the conventional calculation (right plots) 
and the calculation taking into account flavour can be as large as three orders of magnitude 
when Ml = 10^ GeV (see upper left plot) or two orders of mag nitude when Mi = 10^° GeV 
(see lower left plot). 

The reason for this huge enhancement is double. First, in the conventional calculation 
ignoring flavour, when the spectrum has an inverted hierarchy the total CP asymmetry 
goes as Aml^i/ y^Amf^. However, the flavour CP asymmetries go as a/ Am^j^, therefore, 
the individual flavour CP asymmetries can be a factor of 20 bigger than the total CP 
asymmetry. Secondly, the total lepton asymmetry computed ignoring flavour is strongly 
washed-out, since mi > \/ Am^^^. yielding a suppressed baryon asymmetry. On the 
contrary, when there is an approximate texture zero, Aiq, ^ 0, the lepton asymmetry in the 
a-th flavour is only weakly washed-out. These two effects combined are the responsible 
of the huge enhancement of the baryon asymmetry when the spectrum has an inverted 
hierarchy and there is an approximate texture zero in the first row of the neutrino Yukawa 
coupling. 

^^Notice the proximity of both points, which is due to the maximal atmospheric mixing. Furthermore, as 
9i3 approaches zero, the two points coUapse into one, which reflects the fact that in the limit with Oi^ = 
and 023 = 7r/4, there is an exact fj. ^ t symmetry, and imposing A12 = automatically implies A13 = 0. 
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Figure 5: The same as Fig. ([3]), but for a spectrum with inverted hierarchy. 

6 Texture zeros in the two right-handed neutrino 
model 

In this section we would like to study carefully the predictions for the baryon asymmetry in 
the case that there are approximate texture zeros in the neutrino Yukawa coupling. Texture 
zeros commonly arise in model constructions based on the Froggatt-Nielsen mechanism 
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The assignment of different charges under an extra symmetry to particles of different 
generations, translates into a Yukawa coupling with a non-trivial structure in the effective 
theory, once the extra symmetry is spontaneously broken. The assignment of charges could 
be such that the resulting Yukawa coupling could have one or several entries which are very 
small compared to the others. In some instances, the entry for the Yukawa coupling could 
exactly vanish at the high energy scale, although these vanishing entries are usually filled 
when the fields are brought to the basis where the charged lepton Yukawa coupling and 
the right-handed mass matrix are both diagonal. In addition to this, radiative effects can 
also fill these entries. 

In the 2RHN model an exact texture zero in the Yukawa coupling fixes the value of 
tan^ in terms of low energy data. Namely, if \Xii3\ = 0, 

tan^Q^^'' ~ ~^\/ — '"r^ (normal hierarchy) , (68) 

tan^Q^^ ^ —^\/ — "777" (inverted hierarchy). (69) 
V nT'2 

By perturbing the Yukawa coupling around this value of ^q''^ it is possible to lift the 
texture zero, while still reproducing the observed masses and mixing angles. The i?-matrices 
that yield a viable Yukawa coupling with an approximate texture zero are: 

R-[' ^^"C) Vp^^^f ° '"f) 1 (normal hierarchy), 

(70) 

^ , cos^J'^^ ^sin^f^ \ sinOif^ -^cosOU^^ \ ,. ^ , . . , , 

^=1 -sm^r ^cose O)''' [cose}^ esm^f oj (--^^^^ hierarchy), 

(71) 

where pe*'^ parametrizes the departure from the strict texture zero in the ^-parameter space, 
i.e. e = + pe''^. In the case with normal hierarchy the Yukawa couplings explicitly 
read. 



A2. ^y/:^[e-M.^-v/^^pe-(f/:2f/^3-t/^2t^:3)]/t;, (72) 
while in the case with inverted hierarchy, 

^ \fW, [V^^2 (f/*if/^2 - U;,U:,) + ^pe^-Mo.^] /v, 

^^l^,[^M^P-,/w^2pe'-{u:,u;,-u;,u:,)]/v. (73) 

These expressions can be substituted in eqs.(l22l) and eq.(j20l) to derive the CP asymme- 
tries and the washout factors in the a-th flavour, in the case that |Ai/3| ~ 0. Expanding for 
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small values of p and keeping the lowest order terms we obtain for the case with normal 
hierarchy the following CP flavour asymmetries: 



3Mim3 



-Im 





M/3,3 
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f//33 
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and the following washout factors: 
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(75) 



The corresponding formulas for the case with inverted hierarchy are: 
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and 
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(77) 



The predictions for the flavour asymmetries and the washout factors, and accordingly 
for the final baryon asymmetry, will depend on the low energy observables and the "texture 
zero uplifting" parameters p and u. In some cases, the prediction of the baryon asymmetry 
will not depend on p and u, and therefore it could be possible to establish a connection 
between leptogenesis and low energy observables [261 [2H]- Let us analyze separately the 
different possibilites: 



6.1 Texture zero in the (1,1) position 

The possibility of a texture zero in the (1,1) position is perhaps the most interesting one 
from the phenomenological point of view. The postulates that the up and down quark 
matrices are symmetric in the first two generations and that they present a simultaneous 
zero in the (1,1) position, lead to the renown prediction for the Cabibbo angle Ac — 
a/ md/TUg [29]. The success and robustness of this prediction may suggest we apply the 
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same rationale to the leptonic sector, and impose a texture zero in the (1,1) position of the 
neutrino Yukawa matrix, and perhaps also in the charged lepton Yukawa coupling and the 
right-handed neutrino mass matrices. 

The flavour CP asymmetries for the case with normal hierarchy can be straightforwardly 
obtained from eg. (17^. being the result: 

3Mim3 \Me 



2^psin^i3 sin(5 — 0/2 — u) 



Stti;^ 8 \Mee\ 
3Mim3 a/3 ms 

sm 6'i3 



7712 . e- / m2 \ . / e , s 4 



sin 5 + I — sin(5 — 0) H — 1= sin ^13 sin(25 — 
m-i \msj V3 



Sttv^ 8 \M, 



/ \ 2 ^ 



sin (5 + I — sin((5 — 0) ^ sin 6^13 sin(25 — 0) 

?Ti3 V"^3/ V3 



(7^ 

where A^ee — ^36^*^ sin^ ^13 + €^'^1712/ 4: is the (1,1) element of the effective neutrino mass 
matrix (recall that this matrix element is precisely the relevant one for analyses of neutri- 
noless double beta decay). 

On the other hand, the washout parameters are, from eq. ffSll : 

\Mee\ 2 
Kee ^ ^irP ' 

m* 

K,, ^ i^!!^?!^(l + 2v^sin^i3Cos<5), (79) 

with in* ^ 3 X 10~^eV. In view of the present experimental bound sin ^13 < 0.2, it fol- 
lows that lA^eel ^ ■'712/4 and accordingly K^^^rr ^ m3/(2m*) ~ 20. Hence, when there is 
an approximate (1,1) texture zero and p is sufficiently small, the muon and the tau CP 
asymmetries are strongly washed-out, while the electron CP asymmetry is only weakly 
washed-out. 

For the case of an inverted hierarchy of neutrinos, the flavour CP asymmetries read: 

3Mim2 \Mee\ , . 

cosci; sm(0/2j. 



Txv^ ^mim2 2 



e 



3Mim2 3 



~ 87ri;2 



Am^gj sin0 ^ sin^i3 ^2 sin5 + sin(5 — 0) — 3 sin(5 



3Mim2 3 



4 

Am^j„^ a/10 + 6cos0 ' ^3 V V^O + 6cosi 

Am^gj sin sin 6'i3 /2 sin 5 -|- sin((5 — 0) — 3 sin(5 -f- 



87ri;2 8 [Am^^^ VlO + 6cos0 ^3 V VlO + 6cos0 

(80) 

and the washout-parameters, 

l-^eel 2 
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where in this case Aiee — (3mi + e*'^m2)/4. Therefore, it follows that 
Kfj^fj^^rr ^ \/ ^f^ltm/ (2m*), SO the muon and tau asymmetries are strongly washed-out, 
whereas the electron asymmetry is weakly washed-out. 

These formulas can be straightforwardly applied to the recipes presented in subsection 
4.1.3, to go from the flavour asymetries to the baryon asymmetry in the different regimes 
for Ml. We find numerically that when Mi ^10^ GeV it is not possible to reproduce the 
observed baryon asymmetry; in order to reproduce the data it is necessary Mi ^ 10^^ GeV 
so that only the tau Yukawa interactions are in thermal equilibrium. If this is the case, 
leptogenesis only depends on observables that are in principle measurable at low energies, 
since the electron asymmetry is always negligible compared to the muon asymmetry (recall 
that when Mi ^ lO^'^ GeV the relevant quantities to compute the baryon asymmetry are 
^2 = Yee + Y^l^l and K2 = K^e + K^^). Using eqs.( l5T|) . (1571) and taking into account that 
K^^ ^ Kt-t, it follows that when sin ^13 is large, in the case with normal hierarchy Yg oc 
sin^ ^13 sin(2(5— 0), while in the case with inverted hierarchy the relation is more complicated 
(it goes roughly as the term propotional to sin ^13 in eq. flSOj) ). On the other hand, when 
sin ^13 is small, in the case with normal hierarchy the baryon asymmetry is very suppressed 
(it goes as sin^^is), while in the case with inverted hierarchy Yg oc sin0/-\/lO + 6cos(/), 
being in this case the asymmetry suppressed by Am^^j/Am^j^. 



6.2 Texture zero in the (1,2) position 

In the case that there is an approximate texture zero in the (1,2) position of the neutrino 
Yukawa matrix, the flavour CP asymmetries read: 
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(82) 



On the other hand, the washout parameters read: 



Kee — -^(1 — 2A/3sin6'i3 cos5), 
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Therefore, the electron asymmetry and the muon asymmetries are only weakly washed-out, 
while the tau asymmetry is strongly washed-out (when Mi ^ lO^^GeV, I2 would be weakly 
washed-out and Y^^, strongly washed-out). 

We find that when there is an approximate texture zero in the (1,2) position. 
Ml ^ 10^° GeV is necessary in order to reproduce the observed baryon asymmetry. In 
the strict texture zero limit, the baryon asymmetry is dominated by the tau lepton asym- 
metry and therefore there is a well defined connection between leptogenesis and low energy 
observables, Yg oc sin 0. Despite this connection it becomes more diffuse as we depart from 
the texture zero limit, the connection still holds in the region in the vicinity of the (1,2) 
texture zero where the baryon asymmetry is enhanced (see Fig.(jlj) lower left plot). 

We find a similar behaviour when the spectrum presents an inverted hierarchy. In this 
case the fiavour asymmetries are: 
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and the washout factors, 
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with M.^^ — rrii/S + 3/8m2e*'^. Using that lA^^^^I ^ ''^2/2, we find that the electron 
asymmetry is necessarily strongly washed-out, whereas the muon and the tau asymmetries 
are only weakly washed-out (when Mi ^ 10^'^GeV, Y2 would be strongly washed-out and 
Yrr, weakly washed-out). 

For the case with inverted hierarchy, we require again Mi ^ lO^'^GeV to reproduce the 
observed asymmetry (or even larger, when sin ^13 is small). The baryon asymmetry is 
in this case also dominated by the tau asymmetry, except when sin ^13 is very small. In 
the case that the baryon asymmetry is dominated by the tau asymmetry, although there 
exists a connection between leptogenesis and low energy observables, this connection is too 
complicated to be of any practical use. On the other hand, when sin ^13 is very small, there 
is no relation whatsoever, since leptogenesis would depend on the unobservable parameters 
p and uj. 
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6.3 Texture zero in the (1,3) position 



Finally, in the case that the texture zero appears in the (1,3) position, the flavour CP 
asymmetries are: 
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Therefore, in this case, the electron and the tau asymmetries are weakly washed-out, and 
the muon asymmetry, strongly washed-out. On the other hand, for Mi ^ lO^GeV the 
relevant quantity to estimate the washout is K2 = Kee + K^^j^ > 1, so in this regime I2 is 
strongly washed out and Y^t is weakly washed out. 

Similarly to the case of the (1,2) texture zero, this case requires Mi > 10^° GeV to 
reproduce the observations. Furthermore, the baryon asymmetry in the vicinity of the 
texture zero is dominated by the muon asymmetry and hence depends mainly on sin (p. 
This behaviour occurs in particular, in the region where the baryon asymmetry is enhanced 
in Fig.(jlj), lower left plot. 

The case with inverted hierarchy presents some qualitative differences with respect to 
the case with normal hierarchy. When neutrinos have an inverted hierarchy, the flavour 
CP asymmetries read: 
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(89) 



with M.TT — fni/8 + 3/8m2e*'^. As for the case with the (1,2) texture zero, we find that 
the electron asymmetry is necessarily strongly washed-out, whereas the muon and the tau 
asymmetries are only weakly washed-out (also, in the regime where only the tau Yukawa 
interactions are in equilibrium, Y2 would be strongly washed-out and Yrr, weakly washed- 
out). 

It is important to note that all the fiavour asymmetries in eg. (188 p have a suppression 
factor, with different origins. As a consequence, in the limit of the strict texture zero, the 
resulting baryon asymmetry is very small and Mi ^ 5 x lO^^GeV would be necessary to 
accommodate observationl data. However, as we depart from the texture zero limit, we 
find a huge enhancement of the baryon asymmetry, that can allow right-handed neutrino 
masses as low as Mi ~ lO^^GeV, independently of the value of sin ^13. The reason is that 
the tau CP asymmetry can become less suppressed, and at the same time the resulting tau 
lepton asymmetry can be sizable since the tau asymmetry is only weakly washed-out (on 
the contrary, 62 = ^ee + ^^i^J. could be comparable to CrT, but is strongly washed-out). As a 
result, the baryon asymmetry is dominated by the tau asymmetry and hence any connection 
between leptogenesis and low energy observables is lost in this region of enhanced baryon 
asymmetry. 

7 The case of R real 

In Section 5 it was discussed that there are two situations where the differences between 
the computation of the baryon asymmetry taking into account flavour or not are maxi- 
mal, namely when there is an approximate texture zero in the neutrino Yukawa matrix, 
and when the matrix R is real. The case with R real physically corresponds to the class 
of models where CP is an exact symmetry in the right-handed neutrino sector. The rea- 
son for this can be more easily understood working in the basis where the charged lepton 
Yukawa coupling and the right-handed mass matrix are diagonal, so that the neutrino 
Yukawa matrix is the only coupling in the leptonic Lagrangian that violates CP. More 
specifically, the neutrino Yukawa coupling can be written in its singular value decom- 
position, A = Vf[Diag(Ai, A2, A3)Vl. Hence, the CP violation in the right-handed neu- 
trino sector is encoded in the phases in Vr, that can be extracted from diagonalizing the 
combination AA^ = V^Diag(A^, Ag, Ag)!/^. On the other hand, using the parametrization 
of the Yukawa coupling in eq.( l58l) . this same combination of matrices can be written as 
AA^ = M^/"^ RmK^ M^/"^ / v"^ . Comparing the two expressions it is apparent that R is real if 
and only if Vr is real, i.e. when there is no CP violation in the right-handed sector!^ 

In this limit the fiavour CP asymmetries and the baryon asymmetry depend exclusively 

^^Furthermore, it can be checked that there is mixing in Vr if and only if there is mixing in R, and that 
mixing in any 2x2 block in R translates into mixing in the same block of the matrix Vr. 



32 



on the phases of the left-handed sector, that are in turn uniquely determined by the low 
energy phases. Consequently, in this limit the leptogenesis mechanism is tightly connected 
to the low energy phases. This connection is more apparent from the expression of the 
flavour CP asymmetries in the parametrization eq.( l58l) : 

3Mi Y.pJ2j>i3 ^rnfi'm^{m^ - m^)i?i^i?i^Im;7*^;7„^ 

This discussion suggests that the observation of low energy CP violation would consti- 
tute an important hint to the leptogenesis mechanism. In a general case with R complex, 
the low energy phases in the leptonic mixing matrix could stem from the phases in the 
left-handed sector, in the right-handed sector, or in both sectors. In any of the cases, and 
barring unnatural cancellations, a baryon asymmetry is necessarily generated through the 
mechanism of leptogenesis, as long as at least one of the lepton Yukawa interactions is 
in equilibrium (corresponding roughly to M\ < lO^^GeV). This result only follows when 
flavour is correctly taken into account in the Boltzmann equations. In previous analyses of 
leptogenesis ignoring flavour, the observation of low energy CP violation did not automat- 
ically imply the existence of a baryon asymmetry, since the possibility existed that the low 
energy phases could stem exclusively from the left-handed sector and hence be irrelevant 
for leptogenesis. 



8 Conclusions 

Thermal leptogenesis is an attractive and minimal mechanism to make the baryon asym- 
metry of the Universe. The asymmetry is commonly calculated by solving a Boltzmann 
equation for the total lepton asymmetry (one- flavour approximation). In a previous pa- 
per [I3] we studied the impact of lepton flavours (charged lepton Yukawa couplings) on 
the Boltzmann equations (one for each lepton flavour) and discussed the phenomenological 
implications for leptogenesis. 

It may be counter-intuitive that flavour matters in leptogenesis, since Yukawa couplings 
are a small perturbative correction. We have shown that flavour effects are relevant when 
the interaction rates mediated by the the charged Yukawa couplings are faster than the 
typical timescale for leptogenesis. The charged Yukawa rates may be dropped from the 
Boltzmann equations 'provided the latter are written in the flavour basis, where the charged 
Yukawa couplings cannot change the flavour of the asymmetries. This implies that one 
should solve Boltzmann equations for each flavour. 

The final value of the baryon asymmetry depends on the CP asymmetry in each flavour 
a and on the washing out by the lepton number a violating processes. Taking into account 
these flavour dependent washing out factors generically enhances the baryon asymmetry 
with respect to the usual one- flavour approximation, in the limit of strong washout. 



33 



In this paper we have provided analytical approximations for the final baryon asymmetry 
with fiavours accounted for. These depend on the temperature of leptogenesis, and can be 
obtained following the procedure of section 14.1.4^ or of the end of section 14.21 We also 
included CP violation in the AL = 1 scatterings relevant for A'^i production. 

In the two right-handed neutrino (2RHN) model, we have compared our results ob- 
tained with fiavoured Boltzmann equations against the usual one-fiavour approximation, 
to illustrate the big impact that fiavour has on leptogenesis. We have found that there are 
two situations where the differences between the treatment of leptogenesis taking fiavour 
properly into account and previous analyses, that ignored flavour, are maximal. The flrst 
one arises when the neutrino Yukawa coupling present approximate texture zeros in the 
flrst row, so that the CP asjTiimetry in that flavour is only weakly washed-out. As a 
consequence, we have found that thermal leptogenesis in the 2RHN model can produce 
the observed baryon asymmetry for masses of the lightest right-handed neutrino smaller 
than previously believed, namely 10^° GeV for the case with normal hierarchy and 5 x 10^° 
GeV for the case with inverted hierarchy (to be compared with 10^^ GeV and 10^^ GeV, 
respectively, from the conventional computation ignoring flavour). The second situation 
corresponds to the limit in which CP is an exact symmetry in the right-handed neutrino 
sector. In this case, the conventional computation would yield an exactly vanishing baryon 
asymmetry, whereas the computation that takes flavour into account could predict a sizable 
baryon asymmetry. 



Acknowledgements 

It is a pleasure to thank E. Nardi for sharing with us his results and for many useful dis- 
cussions. We also thank L. Covi, G.F. Giudice, A. Pilaftsis and A. Strumia for discussions. 
SD thanks in particular A. Strumia for many thought-provoking and enlightening emails, 
and A. Notari for related results. AR thanks M. Passera for useful discussions. AA and 
FXJM acknowledge the support of the Agence Nationale de la Recherche ANR through 
the project JC05-43009, NEUPAC. 



References 

[1] A.D. Sakharov. JETP Lett., 5:24, 1967. For a review, see A. Riotto and M. Trodden, 
Ann. Rev. Nucl. Part. Sci. 49, 35 (1999) |arXiv:hep-ph/ 9901362|. 

[2] M. Fukugita and T. Yanagida, Phys. Lett. B 174, 45 (1986). 

[3] V.A. Kuzmin, V.A. Rubakov, and M.E. Shaposhnikov. Phys. Lett, B155:36, 1985. 

[4] P. Minkowski, Phys. Lett. B 67 (1977) 421; M. Gell-Mann, P. Ramond and R. Slan- 
sky. Proceedings of the Supergravity Stony Brook Workshop, New York 1979, eds. P. 



34 



Van Nieuwenhuizen and D. Freedman; T. Yanagida, Proceedinds of the Workshop on 
Unified Theories and Baryon Number in the Universe, Tsukuba, Japan 1979, ed.s 
A. Sawada and A. Sugamoto; R. N. Mohapatra, G. Senjanovic, Phys. Rev. Lett. 44 
(1980)912. 

[5] See, e.g. G. F. Giudice, A. Notari, M. Raidal, A. Riotto and A. Strumia, Nucl. Phys. 
B 685, 89 (2004) [arXiv:hep-ph/0310123] . 

[6] W. BuchmuUer, P. P i Bari and M. Plumacher, Annals Phys. 315 (2005) 305 
|arXiv:hep-ph/0401240 . 



[7] A partial list: W. BuchmuUer, P. Di Bari and M. Plumacher, Nucl. Phys. B 



643 (2002) 367 arXiv:hep-ph/0205349]; J. R. Ellis, M. Raidal and T. Yanagida, 



Phys. Lett. B 546 (2002) 228 |arXiv: hep-ph/0206300j ; G. C. Branco, R. Gon- 
zalez Fehpe, F. R. Jo aquim and M. N. Rebelo, Nucl. Phys. B 640 (2002) 202 
|arXiv:hep-ph/0202030|; R. N. Mohapa tra, S. Nasri and H. B. Yu, Phys. Lett. B 
615 (2005) 231 [arXiv:hep-ph/05Q2026]; A. Broncano, M. B. Gavela and E. Jenkins, 
Nucl. Phys. B 672 (2003) 163 [arXiv:hep-ph/0307058l; A. Pilaftsis, Phys. Rev. D 56 
(1997) 5431 |arXiv:hep-ph/9707235|; E. Nezri and J. Orloff, JHEP 0304 (2003) 020 
arXiv:hep-ph/0004227|. S. Davidson and A. Ibarra, Nucl. Phys. B 648, 345 (2003) 



arXiv:hep-ph/0206304|. S. Davidson, JHEP 0303 (2003) 037 arXiv:hep-ph/0302075 



S. T . Petcov, W. Rodejohan n, T. Shindou and Y. Takanishi, Nucl. Phys. B 739 (2006) 
208 |arXiv:hep-ph/0510404|. S. Lavignac, I. Masina and C. A. Savoy, Nucl. Phys. B 633 
(2002) 139 |arXiv:hep-ph/0202086|. E. K. Akhmedov, M. Frigerio and A. Y. Smirnov, 
JHEP 0309 (2003) 021 |arXiv:hep-ph/0305 322]. F. Deppisch, H. Pas , A. Redelbach 
and R. Ruckl, Phys. Rev. D 73 (20'06) 033004 |arXiv:hep-ph/0511062]. 



[8] R. Barbieri, P. Creminelh, A. Strumia and N. Tetradis, Nucl. Phys. B 575 (2000) 61 
[arXiv:hep-ph/9911315 |. 

[9] A. Pilaftsis and T. E. J. Underwood, Phys. Rev. D 72 (2005) 113001 
|arXiv:h ep-ph/050610 7j. A. An isimov, A. Broncano and M. Plumacher, Nucl. Phys. B 
737 (2006) 176 [arXiv:hep-ph/0511248] . 

[10] T. Endoh, T. Morozumi and Z. h. Xiong, Prog. Theor. Phys. Ill, 123 (2004) 
[arXiv:hep-ph/0308276j ; T. Fujihara, S. Kaneko, S. Kang, D. Kimura, T. Morozumi 
and M. Tanimoto, Phys. Rev. D 72, 016006 (2005) |arXiv:hep-ph/0505076j . 

[11] O. Vives, Phys. Rev. D 73 (2006) 073006 |arXiv:hep-ph/0512160| . 

[12] A. Abada and M. Losada, Nucl. Phys. B 673 (2003) 319 |arXiv:hep-ph /0306180]. 
A. Abada, H. Aissaoui and M. Losada, Nucl. Phys. B 728 (2005) 55 
|arXiv:hep-ph/0409343j . 

[13] A. Abada, S. Davidson, F. X. Josse-Michaux, M. Losada and A. Riotto, JCAP 0604 
(2006) 004 |arXiv:hep-ph/0601083| . 



35 



[14] E. Nardi, Y. Nir, E. Roulet and J. Racker, JHEP 0601, 164 (2006) 
|arXiv:hep-ph/0601084| . 

J. A. Harvey and M. S. Turner, Phys. Rev. D 42 (1990) 3344. 

L. Covi, E. Roulet and F. Vissani, Phys. Lett. B 384 (1996) 169 
|arXiv:hep-ph/ 9605319|. 

J. M. Cline, K. Kainulainen and K. A. Olive, Phys. Rev. D 49 (1994) 6394 
| arXiv:hep-ph/9401208] . 

S. Davidson and A. Ibarra, Phys. Lett. B 535 (2002) 25. K. Hamaguchi, H. Murayama 
and T. Yanagida, Phys. Rev. D 65 (2002) 043512 |a rXiv:hep-ph/0109030] . 

J. A. Casas and A. Ibarra, Nucl. Phys. B 618 (2001) 171 |arXiv:hep-ph /0103065]. 

A. Pilaftsis and T. E. J. Underwood, Nucl. Phys. B 692 (2004) 303 
|arXiv:hep-ph/0309342|. 



[15 
[16 

[17; 

[18 

[19 
[20 

[21 

[22 

[23 

[24 

[25 
[26 

[27 
[28 



E. Roulet, L. Covi and F. Vissani, Phys. Lett. B 424, 101 (1998) 
|arXiv:hep-ph/9712468|. W. Buchmuller and M. Plumacher, Phys. Lett. B 431 (1998) 
354 |arXiv:hep-ph/9710460] . 

E. W. Kolb and S. Wolfram, Nucl. Phys. B 172, 224 (1980) [Erratum-ibid. B 195, 
542 (1982)]. 

M. Laine an d M. E. Shaposhnikov, Phys. Rev. D 61 (2000) 117302 
[arXiv:hep-ph799 11473] . 

E. Nardi, Y . Nir, J. Racker and E. Roulet, JHEP 0601 (2006) 068 
|arXiv:hep-ph705 12052] . 

M. Maltoni, T. Schwetz, M. A. Tortola and J. W. F. Valle, New J. Phys. 6 (2004) 122 

A. Ibarra and G. G. Ross, Phys. Lett. B 575 (2003) 279, Phys. Lett. B 591 (2004) 
285 see also Petcov et al., in reference [7j. 

C. D. Froggatt and H. B. Nielsen, Nucl. Phys. B 147 (1979) 277. 

P. H. Frampton, S. L. Glashow and T. Yanagida, Phys. Lett. B 548 (2002) 119. 
M. Raidal and A. Strumia, Phys. Lett. B 553, 72 (2003). S. F. King, Phys. Rev. D 67 
(2003) 113010 [arXiv:hep-ph/0211228|. 



[29] R. Gatto, G. Sartori and M. Tonin, Phys. Lett. B 28 (1968) 128. 



36 



